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Abstract 



In this paper we study generalized classes of volume preserving multidimensional inte- 
grable systems via Nambu-Poisson mechanics. These integrable systems belong to the 
same class of dispersionless KP type equation. Hence they bear a close resemblance 
to the self dual Einstein equation. All these dispersionless KP and dToda type equa- 
r/\ • tions can be studied via twistor geometry, by using the method of Gindikin's pencil of 

two forms. Following this approach we study the twistor construction of our volume 
preserving systems. 

Dedicated to the memory of Dr. B C Guha 



^ ■ 1 Introduction 

In this article we study volume preserving diffeomorphic integrable hierarchy of three 
flows [1]. This is different from the usual two flows cases, and this can be studied via 
Nambu-Poisson geometry. It is already known that a group of volume preserving diffeo- 
morphisms in three dimension plays a key role in an Einstein-Maxwell theory where the 
Weyl curvature is self-dual and the Maxwell tensor is algebraically anti self-dual. Later 
Takasaki [2] explicitly showed how volume preserving diffeomorphisms arises in integrable 
deformations of self-dual gravity. 

Nambu mechanics is a generalization of classical Hamiltonian mechanics, introduced 
by Yoichiro Nambu [3]. At the begining he wanted to formulate a statistical mechanics 
on M 3 , emphasizing that the only feature of Hamiltonian mechanics one should preserve 
is the Liouville theorem. He considered the following equations of motion 

dr _ 

— = Vn(r) A V-y(r), r = (x, y, z) G R, 

where x, y, z are dynamical variables and u, v are two functions of r. Then Liouville 
theorem follows from the identity 

V • (Vw(r) A Vv(r)) = 0. 
Copyright © 2001 by P Guha 
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He further observed from above equation that the evolution of a function on M 3 is given 
by 

df d(f, u, v) 



dt d(x,y, z) ' 

where the right hand side is the Jacobian of the mapping ]R 3 — ► M 3 given by 

(x,y,z) — ► (f,u,v). 

The right hand side can be interpreted as a generalized Poisson bracket. Hence the binary 
operation of Poisson bracket of Hamiltonian mechanics is generalized to n-ary operation 
in Nambu mechanics. Recently Takhtajan [4, 5] has formulated its basic principles in an 
invariant geometrical form similar to that of Hamiltonian mechanics. 

In this paper, we shall use Nambu mechanics to study generalized volume preserving 
diffeomorphic integrable hierarchy. These classes of integrable systems are closely related 
to the self dual Einstein equation, dispersionless KP equations etc. In fact we obtain a 
higher dimensional analogue of all these systems. It turns out that all these systems can 
be written in the following form: 

dnw = o, n {n) a n< n > = o. 

For n = 2 we obtain all the self dual Einstein and dKP type equations. 

Hence we obtain a common structure behind all these integrable system, so there is 
a consistent and coherent way to describe all these systems. Here we unify all these 
classes of integrable systems by Gindikin's pencil or bundle of forms and Riemann-Hilbert 
problem (twistor description) . Gindikin introduced these technique to study the geometry 
of the solution of self dual Einstein equations. Later Takasaki-Takebe [6, 7] applied it to 
dispersionless KP and Toda equations. 

This paper is organized as follows. In Section 2 we present a tacit introduction of 
Nambu-Poisson geometry. In Section 3 we construct our volume preserving integrable 
systems via Nambu-Poisson geometry. If one carefully analyse these set of equations, 
then one must admit that they bear a close resemblance with the volume preserving KP 
equation, so far nobody knows about this equation. It is known that area preserving KP 
hierarchy (= dispersionless KP hierarchy) plays an important role in topological minimal 
models (Landau-Ginzburg description of the A- type minimal models). So we expect 
volume preserving KP hierarchy may play a big role in low dimensional quantum field 
theories. Section 4 is dedicated to the twistor construction of these systems. 

Notations 

1. A4 : Nambu-Poisson manifold. 

2. rj : Nambu-Poisson poly vector (tensor). 

3. J v : Bundle map associated to n. 

4. Xf 1 ...f n _ l : Hamiltonian vector field associated to rj. 

5. L, M, N : Formal Laurent series, Lax operators involving volume preserving sys- 
tems. 



Volume Preserving Multidimensional Integrable Systems 327 

6. C : Lax operator for dispersionless KP hierarchy. 

7. /C : Orlov operator in dispersionless KP. 

8. D, : two form in dispersionless KP. 

9. Q^> : three form in volume preserving integrable systems. 

10. T : curved twistor space. 

11. (P,Q) or (P',Q') : Pair of Darboux coordinates. 

12. e 1 - 7 : one forms or tetrad. 

2 Nambu- Poisson manifolds 

The modern concept of Nambu-Poisson structure was proposed by Takhtajan in 1994 
in order to find an axiomatic formulation for the n-bracket operation. Let Ai denote a 
smooth n-dimensional manifold and C°°(M) the algebra of infinitely differentiable real 
valued functions on Ai . A manifold A4 is called a Nambu-Poisson manifold if there exits 
a M-multi-linear map 

{,..., } : [C°°(M)f n -> C°°(M). 

This is called Nambu-Poisson bracket of order n V /i,/2,--- , iW-i £ C°°(Ai). This 
bracket satisfies 

{/l,...,/n} = (-l) £(ff) *(1)'" ->/«r(n)}. (2- 1 ) 

{/l/2, /3, • • • 1 /n+l} = /l{/2, /3) • • • j fn+l} + {/l, /3, • • ■ j /n+l}/2, (2.2) 

and 

{{/lj--- , fn-l, fn}, fn+l, ■ ■ ■ , /2n-l} 

+ {/n,{/l,--- , /n-li/n+l},/n+2, • • • j/2n-l} + - " 

+ {/n, • • • j jln-2i {/l, • • • j /n-1) /2n-l}} 

= {/l> • • • j /n-1, {/n, • • • , /2n-l}}, (2.3) 

where a € S n -- the symmetric group of n elements -- and e(<r) is its parity. Equa- 
tions ( |2.lD and fl2.2| ) are the standard skew-symmetry and derivation properties found for 
the ordinary (n = 2) Poisson bracket, whereas ( |2.3j ) is a generalization of the Jacobi iden- 
tity and was called in [2] the fundamental identity. When n = 3 this fundamental identity 
reduces to 

{{/l,/2,/3},/4,/5} + {/3, {/l,/2,/4},/5} 

+ {/3,/4, {/l,/2, /5}} = {/l,/2, {/3, /4, /5}}- 

2.1 Hamiltonian geometry 

The Nambu-Poisson bracket is geometrically realized by the Nambu-Poisson polyvector 
77 e r(A"TM), a section of ATM, such that 

{/l,--- ,/n} =T](dfi,... ,df n ), 
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in local coordinates (xi, . . . , x n ) this is given by 
, , d d 



OXi^ OXi n 

where summation over repeated indices are assumed. 

The structure is regular if 77 7^ 0. The canonical example of a Nambu-Poisson structure 
of order n > 2 is the one induced by a volume form on an oriented manifold 

{/i,--- ,/n} = v(dfi,--- ,df n ). 

The Nambu-Poisson polyvector defines a bundle map 

J„ : T*A4 x • • • T*M — ► TM 

' * ' 

n— 1 

given by 

(P,J v (ai,... ,a n _i)) = r/(ai,... ,a n -i,/3), 

where /3,a; <E O^M). 

The Hamiltonian vector field is defined by 

X fl,...,f n -1 = Jri{dfl, ■ ■ ■ ,df n -l)- 

This is a Hamiltonian vector field of (n — 1) functions. Thus, we can express the funda- 
mental identity as 

Lx fi~f„-M d 9i ■ ■ ■ dg n ) = 0, 

where L is the Lie derivative. This shows that the Hamiltonian vector fields are the 
infinitesimal automorphism of the Nambu-Poisson tensor. 

The polyvector field rj € T{f\ n TM) defines an n-ary Poisson bracket if and only if 
either n is even and the Schouten bracket satisfies [77, rj\ = 0, or n is odd if rj satisfies the 
conditions (see for detail [81) 



n , 

iafa)Ai0(r7) = O, Va,[3£T*M, ^(idxjV) A lL 

.7=1 ^ 



_g_»/) =0. 



In this framework we can introduce the bihamiltonian structure also. A bihamiltonian 
system is prescribed by specifying two sets of Hamiltonian function H = Hi, . . . , H n and 
h = hi,... ,h n , where H { , hi G C°°(M), 

X = J°(dhi,. . . , dh n ) = JJjidHu ... , dH n ), 

where JIs are the bundle maps. 
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One more interesting algebraic structure pops up in the Nambu-Poisson geometry - 
this is called Leibniz algebra. The Nambu-Poisson polyvector rj induces a homomorphism 
of vector bundles 

J* : A fe (T*M) — ► A n " fc (TM), 

such that J~(f3) = ipr} for (3 € f\ k T*M. All these higher covectors form Leibniz algebra. 
It was shown by Ibanez et. al. [9] that the bracket of (n — 1) forms on A4 is the M-bilinear 
operations 

[vlLdtata = C 0O (A n - 1 T*>f) x C°°(A n - 1 T*7W) — ► C°°(A n - 1 T*M), 
given by 

[a,[3)=L Ja p + {-l) n {i{da)r ] )f3, 

for all at,P € C°°(A n - 1 T*7W). 

Any Leibniz algebra (A, o) satisfies 

fo(goh) = (fog)oh + go(foh), f,g,h£C°°(M). 

A morphism in the Nambu-Poisson category is a map 

<f> : M — > N 
between Nambu-Poisson manifolds preserving Nambu-Poisson brackets: 
{/i, • • • , Iu}n ° = {/i o 0, • • • , /„ o </>} M . 

The fundamental difference between the Nambu-Poisson bracket and the classical Pois- 
son case is that for n > 3 the Nambu-Poisson polyvector r] is decomposable, i.e. it has 
rank n at points where it does not vanish. 

The Nambu-Poisson polyvectors were charcterized by Takhtajan. The polyvector 
field rj is a Nambu-Poisson tensor if and only if the natural component of r] satisfy certain 
conditions. It was stated in [4] that the fundamental identity fl2.3|) is equivalent to the 
following algebraic and differential constraints on the Nambu-Poisson tensor rj: 

(A) Sij + P(S) i3 = 0, 

for all multi-indices i = {ii, . . . ,i n } and j = {ji, ... ,j n } from the set {1, . . . , iV}, where 

&ij = Vh—inVjl—jn "i Vjnilis—inVji—j„-ii2 + ' ' ' 

' Vjn.i2—in-lilVjl»-jn-lin ~ Vjn.i2—in.Vjl—jn-lh'l 

(B) P is the permutation operator which interchanges the indices i\ and j\ of 2n- 
tensor 5, and 

/ . I *2-«n Q "T 'Ijnh3-ln q 1 I" >lj n l2...in-ll a I 

T 1 \ * « " / 



Vjljl-jn-ll' 



dxi 
1=1 ' 

for all z 2) ... ,i n ,ji,... ,j n = 1, . . . ,iV. 
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It was conjectured in [2] that the equation Sij = is equivalent to the condition that 
n-tensor r\ is decomposable, recently this has been proved by Alekseevsky and author [10], 
Since then several other proofs appeared this surprising result (see [11, 12] and references 
therein). 

The algebraic identity is the special feature for the Nambu-Poisson structure. It ap- 
pears solely in the Nambu-Poisson geometry — there is no counterpart of this structure 
in the usual Poisson geometry. 

The following theorem describes the local structure of the Nambu-Poisson structure of 
order n with n > 3. 

Theorem 2.1 (Local Triviality Theorem). Let r/ be a Nambu-Poisson polyvector of 
order n > 3. Near any point at which rj does not vanish there are local coordinates 
xi, . . . ,x n such that 

d d d 

V = t— A -— A • • • A 



dx\ 8x2 dx 



n 



Let {hi, . . . , h n -i} be the set of Hamiltonian functions, and Xf ll ...h„_ 1 be corresponding 
the Hamiltonian vector field. Then the equation for the integral curves of Xhi—hn-i is 

x = J v (dhi,... ,dh n -i), 

Hamilton's equations of motion. 

2.2 Nambu-Poisson dynamics 

It is also shown in [2] that Nambu dynamics on a Nambu-Poisson phase space involves 
n — 1 so-called Nambu-Hamiltonians Hi,... ,H n ^i £ C°°(Ai) and is governed by the 
following equations of motion 

f t ={f,Hi,...,H n _i}, V/GC°°(X). 

A solution to the Nambu-Hamilton equations of motion produces an evolution operator 
Ut which by virtue of the fundamental identity preserves the Nambu bracket structure on 

C°°{M). 

Definition 2.1. / £ C°°(M) is called an integral of motion for the system if it satisfies 
{/) Hi,H 2 , . . . , i? n -i} = 0. 

Example. Let us illustrate how Nambu-Poisson mechanics works in practise. The exam- 
ple is the motion of a rigid body with a torque about the major axis introduced by Bloch 
and Marsden [13]. 

Euler's equation for the rigid body with a single torque u about its major axis is given 
by 

rhi = aim2'm'3, w-2 = a2m\raj„ 1713 = a^mim2 + u, 

where u = —kmirci2 is the feedback, ai = j — jr, 0,2 = j — 7- and 03 = j- — j-. We 
assume I\ < I2 < I3. 
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These equations can be easily recast into generalized Nambu-Hamiltonian equations of 
motion 

drrii 



{Hi,H 2 ,rrii}, 



dt 

where the right hand side is given by 

(tt tt i d(H 1 ,H 2 ,m i ) 

{Hi,H 2 ,mi} ■= ^7 r- 

o(mi,m 2 ,m3) 

These equations involve two Hamiltonians and these are 

Hi = - [a 2 m l - aim 2 ) , H 2 = - I m 1 - m 3 I . 

When a\ = a 2 = 03 = 1 and 113 = 0, these set of equations reduce to a famous Euler 
equation or Nahm's equation 

~~Tr = e ijk[l ji J-k\i % i3i<Z = 1,2,3, 

where T^s are SU{2) generators. 

It is well known that the system of TV vortices can be described by the following system 
of differential equations [14] 

N 

— 1 



where z n = x n + iy n are complex coordinate of the centre of n-th vortex. 

The equation of motion of a system of three vortices can be put in the following form [15] 

Mi = e i j k Y i M i {M j - M k ), i,j, k = 1, 2, 3 

where Mi = \z 2 — 23 [ 2 , M 2 = \z% — z\\ 2 and M3 = \z\ — z 2 \ 2 . Again, these equations 
can be described by Nambu-Poisson geometry. In fact, there is a large scope to apply 
Nambu-Poisson formulation in fluid mechanics [16]. 

3 Volume preserving integrable systems 

In this section we shall follow the approach of Takasaki-Takebe's [6, 7] method of area 
preserving diffeomorphic (or sDiff(2)) KP equation. In fact, they adopted their method 
from self dual vacuum Einstein equation theory. In the case of self dual vacuum Einstein 
equation, and hyperKahler geometry also area preserving diffeomorphism appear, where 
the spectral variable is merely parameter. But in case of sDiff(2) K.P. equation the situa- 
tion is different, where Takasaki-Takebe showed that one has to treat A as a true variable 
and it enters into the definition of the Poisson bracket. Our construction is closely related 
to Takasaki-Takebe construction, only the Kahler like two form and the associated "Dar- 
boux coordinates" is replaced by volume form and the Poisson bracket is replaced by its 
higher order Poisson bracket called Nambu bracket. 

Suppose we consider L = L(X,p,q), M = M(X,p,q) and N = N(X,p,q) are some 
Laurent series in A with coefficients are functions of p and q. 
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Definition 3.1. The volume preserving integrable hierarchy is defined by 

dt n 
8M 



., {Bi n ,B 2n ,L}, (3.1) 



■ {B ln ,B 2n ,M}, (3.2) 

{B ln ,B 2n ,N}, (3.3) 

and 



dtn 

dN 

dt n 



{L,M,N} = 1, (3.4) 

where B\ n := (L n ) n >o and B 2n '■= (M n ) n >Q. The first three equations are hierarchy 
equations, and these generates the three flows of the system and the last one shows the 
volume preservation condition. The Nambu-Poisson bracket in 3D "phase space" (X,p,q) 
is given by 

dA I dB dC dB dC\ 
{A(X,p,q),B(X,p,q),C{X,p,q)} = — [-q—q ~Q~~d~) + cyclic terms - 

Let us now compare our case with the area preserving KP hierarchy. The sdiff(2) KP 
hierarchy is given by 

^- = {B n ,C}, ^ = {B n ,lC}, {£,£} = 1, 

where £ is a Laurent series in an indeterminant A of the form 

oc 
£ = X + Y,U n+ l(t)X- n , 
n=\ 

B n = (jC)>o- The function /C is called Orlov function and it is defined by 

oo oo 

n=l i=l 

where t\ = x. 

Remark 3.1. Our hierarchy has a structure of volume preserving KP hierarchy and in- 
stead of one Orlov function /C, we need two Orlov functions M and N. 

Any two equations of the hierarchy commute 

dt m d tl =d tl d tm , 

hence we obtain 
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Proposition 3.1. The Lax equation for L, M, N are equivalent to the following equa- 
tions: 

9Bm n 1 , f„ dB 2 



dt m 



,B 2n > + S Bi n , r, \ + \Bi,B 2m > 



+ \B lm ,H 2 } = 0, 
L vz n ) i ui n ) I 

where 

Hi = {Bi n , B 2n , Bim}, H 2 = {Bin, B 2n , B 2m }. 

Proof. Result follows from the compatibility conditions of hierarchy equations and the 
fundamental identity. ■ 

Remark 3.2. In the case of sDiff(2) hierarchy the above expression boils down to zero 
curvature equation. 

Let fty> be a three form given by 

Definition 3.2. 

oo oo 

^(3) . = ^ dBin A dB 2n A dt n = dX A dp A dq + ^ dB Vn A dB 2n A dt n . (3.5) 

n=l n=2 

From the definition it is clear 0, is closed 3 form. In fact sDiff(3) structure is clearly 
exhibited from this structure and the theory is integrable in the sense of nonlinear graviton 
construction [14]. This is a generalization of nonlinear graviton construction. 

Theorem 3.1. The volume preserving hierarchy is equivalent to the exterior differential 
equation 

ft(3) = dL A dM A dN. (3.6) 

Proof. We have seen that £1 can be written in two ways. Expanding both sides of the 
exterior differential equation as linear combinations oid\f\dpf\dq, d\AdpAdt n , d\AdqAdt n 
and dp A dq A dt n . 

When we pick up the coefficients of dX A dp A dq, we obtain the volume preserving 
condition 

{L,M,N} = 1. 

When we equate the other coefficients, viz. dXAdpA dt n , dXAdqA dt n and dpAdqA dt n 
we obtain the following identities: 

d(Bm,B 2n ) _ d(L,M,N) d(Bi n ,B 2n ) _ d(L,M,N) 

d(X,p) d(X,p,t n ) ' d(X,q) d(X,q,t n ) 

and 

d{Bm,B 2n ) _ d(L,M,N) 
d(p,q) d(p,q,t n ) 
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respectively. We multiply the above three equations by |£, Sj and |j respectively. If we 
add first and third equations and substract the second one from them, then after using 
the volume preserving identity we obtain 

-Q— = {Bin, Bin, L}. 

Similarly we can obtain the other equations also, in that case we will multiply the equations 

m dM QT1 j dM_ 
dq ' dp dUU ax 



(O), Q and (O) by m, ™ and m respectively. 



Equation (^) and (|3 . 6|) , we obtain 

d f Af dL AdN + ^T B ln dB 2n A din J = 0. 
This implies the existence of one form Q such that 

oo 

dQ = Md(LdN) + Y, B ln d(B 2n dt n ). 



n=l 



This is an analogue of "Krichever potential" in the volume preserving case. Hence we 
can say from (|3.4[) 



M dQ I B 9Q I 

8{LdN) lB2n ' tn fixed ' ln 9(S 2n dt n ) lL,Ar ' B2m ' im(m ^ n)fixed ' 

4 Application to mulitidimesional integrable systems 
and Riemann— Hilbert problem 

We already stated that our situation is quite similar to nonlinear graviton construction 
of Penrose [17] for the self dual Einstein equation. This is a generalization of nonlinear 
graviton constructions. 

To the geometer self dual gravity is nothing but Ricci flat Kahler geometry and it is 
characterized by the underlying symmetry groups sDiff(2), this are called area preserving 
diffeomorphism group on surfaces. These are the natural generalization of the groups 
Diff(S' 1 ), diffeomorphism of circle. 

It is well known how the area preserving diffeomorphism group appears in the self dual 
gravity equation. Let us give a very rapid description of this. 

Let us start from a complexified metric of the following form 

j 2 j j. I e e \ 11 22 12 21 

ds = det ( 21 2 2 I = e e — e e 

where e 13 are independent one forms. Ricci flatness condition boils down to the closed- 
ness of 

dn kl = o (4.i) 
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of the exterior 2-forms 

n kl = ^J ije ik Ae jl , 

where J is the normalized symplectic form 
1 



1 • -1 

Then above system of two forms can be recast to 
fi(A) = i J zj (e il + e l2 X) A (e jl + e j2 A) . 

This satisfies 

0(A) A Q(X) = 0, dfi(A) = 0, (4.2) 

where (i stands for total differentiation. These suggest us to introduce a pair of Darboux 
coordinates 

n{X) =dPAdQ 

and these are the sections of the twistor fibration 

tt:T — >CP\ 

where T is the curved twistor space introduced by Penrose. Basically each fibre is endowed 
with a symplectic form and as the base point moves this also deforms and here comes the 
area preservation. 

Two pairs of Darboux coordiantes are related by 

/(A, P(A), Q(A)) = P', g(X, P(A), Q(A)) = Q' 

and / and g satisfy {/, g) = 1. The pair (/, g) is called twistor data. Locally / and g 
(after twisting with A) yield patching function. Ricci flat Kahler metric is locally encoded 
in this data. This set up is nothing but the Riemann-Hilbert problem in area preserving 
diffeomorphism case. 

The novelty of this approach is that this twistor construction will work in the higher 
dimensions too, when there is no twistor projection. The most important example is 
the electro-vacuum equation, volume preserving diffeomorphism groups in three dimen- 
sion play a vital role here. This model was first introduced by Flaherty [18] and later 
Takasaki [2] showed how this works explicitly. 

4.1 Gindikin's bundle of forms 

We already stated that anti-self dual vacuum equations govern the behaviour of complex 
4-metrics of signature (+,+,—,—) whose Ricci curvature is zero and whose Weyl curva- 
ture is self dual. These two curvatures are independent of change of coordinates, so in 
one particular of the equations these metric becomes autometically Kahler and can be 
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expressed in terms of a single scalar function f2, the Kahler potential. Then curvarure 
conditions will lead you to 1st Plebenski's Heavenly equation 

d 2 n d 2 n d 2 n d 2 n 

= 1, 

dxdx dydy dxdy dydx 

or using the Poisson bracket (w.r.t x and y) 

and the correspoding anti-self-dual Ricci flat metric is 

d 2 n ... 

ff(^) = gxid ~ j dx t dx 3 , x l = x,y, x 3 =x,y. 

This system is completely integrable and it is an example of a multidimensional inte- 
grable system. 

Let Q, be the 2-form, given by 



+ X(Cl xx dx A dx + £l X ydx A dy + & yx dy A dx + VL y ydy A dy) + X 2 dx A dy. 



fl(x, x, y, y) = dx A dy 
Lemma 4.1. 

i) dn = o, 

2) fiAO = 0. 
Proof. Since S7 satisfies Plebenski's Heavenly 1st equation, hence 2) is true. ■ 

A number of multidimensional integrable systems can be written in terms of a 2-form 
f2 which satisfies the equations ( |4.l[) . 

Example. The dispersionless KP hierarchy has a Lax representation with respect to a 
series of independent ("time") variables t = (ti, £2, • • • ) 

dT n ={Bn > C} > 

where B n := (£ n )>o, n= 1,2,..., £ is a Laurent series in an indeterminant A of the form 
00 

n=l 

{ , } is a Poisson bracket in 2D phase space with respect to (A, x). 
Let us consider 

00 
Q, = d\ A dx + ^^ dB n A dt n 

n=2 

then f2 A $7 is equivalent to the zero curvature condition 

oB„ dB rn , „ „ , 
^--^ + {B n ,B m } = 0. 

Morn VLn 

This is an alternative form of dispersionless KP hierarchy. 

All these systems are related to area preserving diffeomorphism group sDiff(2). 
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In this paper we are presenting an analogous picture for multidimensional integrable 
systems related to volume preserving diffeomorphism group. 

Proposition 4.1. The 3- form 

oo 
Q( 3 ) = d\ A dp A dq + ^ dB ln A dB 2n A dt n , 

n=2 

satisfies 

dn {3) = 0, fl (3) A ft (3) = 0. 



Proof. This is equivalent to Propostion 3.1, hence it is satisfied. 



The Gindikin's method [19] of pencil of 2-forms is the most effective way to study these 
systems. Consider the following system of 1st order equation depending on a parameter 

t = (ti,t 2 )gC 2 



e\r) = e^r k 1 +--- + e lk rt 



e 2( r ) = e 2i r fc + ... + e 2fe r . 



2 ! 



e a (r) = e 2,1 7f + • ■ ■ + e 2I M, 



where e ij are 1-forms. Let Q, k (r) be the bundle of 2-forms 

Q k ( T ) = e i( r ) a e 2 (r) + • • • + e^V) A e 2 '(r) 
satisfying the conditions 

(n k ) +1 = o, (n k ) + o, do fc = o. 

The bundle of forms actually encodes the integrability of the original system. In the special 
case I = 1, k = 2, we recover the Ricci flat metric 

5 = e ll e 22_ e 12 e 21^ 

4.2 Higher dimensional analogue of Gindikin's Pencil 

Let us consider the following system of 1st order equations depending on parameter r = 
(ti,t 2 ,t 3 ) GC 3 , 

e x (r) = e n Ti + e 12 r 2 + e 13 T 3 , 
e 2 (r) = e 21 n + e 22 r 2 + e 23 T 3 , 
e 3 (r) = e 31 ri + e 32 r 2 + e 33 T 3 . 

Just like the previous situation these bundle of forms then encodes the integrability of 
the original system. The metric defined here is the higher order analogue of Ricci metric 
in two form case. 

This metric is given by, 

g = e H e 22 e 33 _ e ll e 32 e 23 + ^31^,23 _ ^21^3 + ^21^,32 _ ^31^ 
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Remark 4.1. Since 3-form Q^ 3 ' satisfies 

dn {3) = o, n (3) a n {3) = o, 

so it denotes the volume preserving multidimensional integrable systems. These systems 
can be described by Gindikin's bundle of multi forms, higher dimensional analogue of 
nonlinear graviton. 

4.3 Twistor description of volume preserving multidimensional 
integrable systems 

The natural question would be to find out the analogous Riemann-Hilbert problem in 
the volume preserving case. Our situation is very similar to electro-vacuum equation. 
Let us consider two sets of solutions of hierarchy (L, M, N) and (L, M, N) with different 
analysity. Then there exist an invertible functional relation between these two sets of 
functions such that it satisfies 

L = h(L,M,N), M = f 2 (L,M,N), N = f 3 (L,M,N), 

where /i = fi{X,p,q), h = f2(\p,q) and f 3 = fz{\,p,q) are arbitrary holomorphic 
functions defined in a neighbourhood of A = oo except at A = oo. 
We assume /i,/2,/3 satisfy the canonical Nambu Poisson relation 

{/l,/2,/ 3 } = l- 

This is a kind of Riemann-Hilbert problem related to three dimensional diffeomor- 
phisms. In this case sDiff(3) symmetries is clear, in fact sDiff(3) group acts on (/i, f 2 , fs), 
we can lift this action on (L, M, N) and (L, M, N) via Riemann-Hilbert fatorization. 

4.4 Application to hydrodynamic type systems 

There are certain kind of integrable systems, called hydrodynamic type, naturally arise in 
gas dynamics, hydrodynamics, chemical kinematics and may other situations, these are 
given by 

u\ = v l j (u)u J x , 

where vUu) is an arbitrary N x N matrix function of 

u = (u 1 , . . . ,u N ) , u l = u l (x,t), i = l,...,N. 

The Hamiltonian systems of hydrodynamic type systems considered above have the 
form 

ul = {u\H}, 

where H = j h(u)dx, is a functional of hydrodynamic type. The Poisson bracket of these 
systems has the form 

{u i (x),u j (y)} = g ij (u)S x (x - y) + ^(u)u k J(x - y), 
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called Dubrovin-Novikov type Poisson bracket. Dubrovin-Novikov showed if the metric is 
non degenerate i.e. det [g 4J ] 7^ then the above bracket yields a Poisson bracket provided 
g lJ (u) is a metric of zero Riemannian Curvature. 

A large class of them can be described by the Lax form, 

ip x = zAip, ip t = zBifj. 

From the compatibility condition we obtain, 

A t = B x , AB = BA. 

These set of equations can be easily recasted to 

dn = 0, fiAO = 0, 

where $1 is a closed one form. 

Hence has also a twistorial description. The upshot of this section is that Gindikin 
pencil of forms can be applied to large number of classes integrable systems, including the 
volume preserving integrable systems we propose here. 

5 Conclusion 

In this paper we have shown that there are large class of integrable systems can be ob- 
tained from Nambu-Poisson mechanics. They belong to the same family of self dual 
Einstein or dispersionless KP type equations. In some sense these integrable systems are 
higher dimensional generalization of self dual Einstein equation. Hence these systems are 
describable via Gindikin's bundle of forms, or twistor method. Our integrable systems are 
volume preserving. We hope that eventually such construction will find their use in the 
construction of volume preserving KP and Toda type equation. 

We have also discussed in this paper that Nambu-Poisson manifold is an useful tool to 
study volume preserving integrable systems. Recently in membrane theory physicists [20] 
have found M-algebra from M-brane, these are related to Nambu-Poisson mechanics. 

There are certain problems we have not discussed in this paper, viz. the quantization 
of these volume preserving generalized multidimensional integrable systems. Presumably, 
the method of star product quantization [21, 22, 23] would be the best way quantize these 
systems and instead of binary star product we need triple star product [24] . 
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